Thermally assisted magnetization reversal in the presence of a spin-transfer torque 
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We propose a generalized stochastic Landau-Lifshitz equation and its corresponding Fokker-Planck 
equation for the magnetization dynamics in the presence of spin transfer torques. Since the spin 
transfer torque can pump a magnetic energy into the magnetic system, the equilibrium temperature 
of the magnetic system is ill-defined. We introduce an effective temperature based on a stationary 
solution of the Fokker-Planck equation. In the limit of high energy barriers, the law of thermal 
agitation is derived. We find that the Neel-Brown relaxation formula remains valid as long as we 
replace the temperature by an effective one that is linearly dependent of the spin torque. We carry 
out the numerical integration of the stochastic Landau-Lifshitz equation to support our theory. Our 
results agree with existing experimental data. 
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I. INTRODUCTION 

Thermally assisted magnetization reversal has been 
the subject of intensive theoretical and experimental 
study for many decades. Aside from the relevance of 
the subject to emerging magnetic technology such as 
heat-assisted magnetic recording fill and thermal stabil- 
ity of magnetic random access memory |^], fundamen- 
tal physics of magnetization reversal process driven by 
white-noise is very rich. Classical transition-rate theory 
of Kramer 3] has supplied a framework in understand- 
ing thermal activation of a single domain magnetic ele- 
ment 01 . Namely, the thermal switching probability P{t) 
can be described by the Neel-Brown (NB) relaxation- 
time formula, P{t) = 1 — exp(— i/r), where the relax- 
ation time is r = fQ^exp{Eb/kBT), /o is an attempt 
frequency, Ei, is the energy barrier, and T is the temper- 
ature. For a multi-domain structure, the energy surface 
becomes extremely complicated and identifying energy 
barriers are numerically non-trivial. Nevertheless, with 
recent development of micromagnetic modeling, one can 
understand thermal reversal reasonably well for a not- 
too-complicated structure (H. 

An implicit and yet critical assumption in the NB 
theory is that magnetization dynamics is governed by 
a torque from an effective magnetic field Heff = 
— Vm£^(M), where E{M.) is the total magnetic energy, 
i.e., the effective field is derivable from the derivative 
of an energy function with respect to the magnetization 
vector. Therefore, an energy barrier is well-defined in 
the NB relaxation formula. If the torque is not derivable 
from an energy function, one would expect breakdown of 
the NB formula. Recently, a new class of torques, called 
spin transfer torque (STTL has been proposed |EI3 and 
verified experimentally [^Q- STT is derived from a spin 
polarized current in magnetic multilayers. For a spin 
valve structure, STT is written as 



r, = I^M X (M X Mp) 



(1) 



where aj represents the strength of STT; it is propor- 



tional to the current density. 7 is gyromagnetic ratio, 
Mp is a unit vector representing the direction of the 
magnetization of the pinned layer, M is the magneti- 
zation vector of the free layer and Mg = |M| is the 
saturation magnetization. If we define an effective field 
H;^ = {aj/Ms)M X Mp from STT, it is evidently that 
H J can not be written as a total derivative of a function 
with respect to the magnetization vector, i.e., there is no 
well-defined energy associated with the field Hj. 

Recent experiments 0, 0, ^2 the thermal effect of 
the spin torque had also indicated that the thermally as- 
sisted magnetization reversal can not be simply described 
by the NB formula. Urazhdin et al found that the 
activation energy strongly depends on the magnitude as 
well as the direction of the current. To capture the gloss 
features of the observed experiments, they have to in- 
troduce an effective temperature unrelated to the true 
temperature in the NB formula. Their proposed effective 
temperature was then interpreted via a possible magnetic 
heating and magnetic excitations from the spin transfer 
torque. The current directional dependence of the effec- 
tive temperature indicated that the heating is not of the 
ordinary current-induced Joule heating. However, no at- 
tempt has been made to mathematically link the effective 
temperature with the spin transfer torque of Eq. (1). 

The problem of thermally assisted escape process 
driven by a non-gradient driven force, not derivable from 
a potential, is an unresolved outstanding problem in sta- 
tistical physics. While there are already some efforts to 
formulate the escape time in this case, the general conclu- 
sion is that the law of escape time lacks universality and a 
variety of scaling relations exist ^3 ■ The standard treat- 
ment of the thermal escape problem in the presence of 
a non-gradient field would start from the Fokker-Planck 
(FP) equation and one numerically solves for the prob- 
ability distribution H^. This procedure involves proper 
averaging over the possible escape trajectories. Let us 
consider the total work done by the conservative torque 
and the non-conservative STT along an arbitrary trajec- 
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tory, 



5W 



aj 



(M X e^) • dM 



(2) 



Mo 



where we have assumed that the magnetization vector 
starts at an initial equihbrium point Mq = M^e^ and 
reaches to an energy saddle point M^, and we have de- 
fined the energy barrier from the conservative torque 
Eb — iJ(M/) — i?(Mo). One immediately reahzes that 
6W defined above depends on the escape trajectory. In 
the absence of STT, one relies on the assumption that the 
most probable path of the escape is through a minimum 
energy barrier, i.e., My would be an energy saddle point. 
In the presence of STT, such assumption breaks down in 
general. In the present paper we do not intend to address 
the general problem of the thermal escape in an open sys- 
tem, rather we want to answer a very focused question: 
to what extent one can formulate the thermal agitation in 
terms of a simplified Neel-Brown activation process and 
to what accuracy one can analyze the relevant experi- 
mental data through a simple effective formula we will 
develop in the later sections. The paper is organized as 
follows. In Sec. II, we propose the stochastic Landau- 
Lifshitz equation and its corresponding Fokker-Planck 
equation in the presence of the current. In Sec. Ill, we 
introduce a stationary solution for the probability den- 
sity of magnetization by identifying an effective barrier 
or an effective temperature associated with spin torques. 
In Sec. IV, we present a numerical calculation to demon- 
strate the validity of our theory in several realistic cases. 
Finally, we compare our theory with existing experiments 
and summarize our theory. 



II. STOCHASTIC LANDAU-LIFSHITZ 
EQUATION IN THE PRESENCE OF CURRENTS 

Let us first explicitly propose the following generalized 
stochastic Landau-Lifshitz (LL) equation that describes 
dynamics of the magnetization vector subject to a STT 
at finite temperatures 



^ --7Mx(He//- 



7a 



h,)-^Mx [Mx (He//+hO]+r„ 

(3) 

where a is the damping constant, Heff is the effective 
magnetic field including the external field, the anisotropy 
field, the exchange field, and the demagnetization field, 
and hj. is a fluctuating field with a Gaussian stochastic 
process whose statistical properties are defined as 

< hl.{t) >= 0; < h\.{t)hi{t') 2D5,j5{t - t') (4) 

where i and j are Cartesian indices, D represents the 
strength of the thermal fluctuations whose value will be 
determined later. <> denotes an average taken over 



all realization of the fluctuating field. In the absence 
of the spin torque, the above equation is the standard 
stochastic LL equation. Note that we have conveniently 
dropped the customary renormalized gyromagnetic ratio 
7/(1 + a^) when compared with the standard Landau- 
Lifshitz- Gilbert equation. The critical assumption of our 
proposed stochastic LL equation, Eq. (3), is that the spin 
torque does not contain a fluctuating field h^. The justi- 
fication for this choice is that the spin torque comes from 
the conduction electrons whose transport properties are 
less affected by thermal fluctuations since the Fermi level 
is much higher than the thermal energy. Therefore, the 
thermal fluctuation would not appear to affect aj which 
represents the strength of the spin torque. We believe 
that our proposed stochastic LL equation captures the 
main random processes induced by thermal fluctuation. 
Nevertheless, one could, in principle, have introduced a 
second random field or torques to take into account the 
fiuctuation of the spin torque. In our proposed LL equa- 
tion, the thermal effect on the spin torque is only encoded 
in the dependence of the magnetization vector. 

To establish the thermal properties from the above 
stochastic equation, one must take a proper thermal av- 
erage. Fortunately, much of theoretical work on the 
stochastic LL in the absence of the spin torque had been 
carried out 0, Q| . Here we will follow and generalize the 
procedure pioneered by Brown We define -P(M, t) as 
a non-equilibrium probability density for magnetic ori- 
entation vectors associated with the stochastic Eq. (3). 
The rate equation for P(M,t) is 



dP 
'dt 



V • J - AV^P = 



(5) 



where the probability current density J = PdlsA/dt, V 
is a short notation for the gradient operator on the mag- 
netization vector Vm, and A is the diffusion constant 
whose value is determined by fluctuation-dissipation the- 
orem. In the present case, A is related to by A = 
(l/2)£)7o(l + a^). The above rate equation, Eq. (5), is 
a simple statement for the angular momentum conser- 
vation: the change of probability density in an enclosed 
small volume (first term) has to be balanced by the net 
probability in-flowing flux (second term) plus the prob- 
ability density loss via spin diffusion (third term). By 
inserting Eq. (3) into Eq. (5), after a straightforward but 
rather tedious algebra manipulation [T5l |. the resulting 
equation is the Fokker-Planck equation 



^ _V • I [-7M X He// + r, - X (M X He//) 

-^7^(1 + a^)DM X (M X V)] P} 

(6) 

In the absence of the spin torque, i.e., Tg = 0, the thermal 
equilibrium distribution density P demands to take the 
form of the Boltzmann distribution function, i.e., P{aj = 
0, T) cx exp(— iJ/fc^r) where T is the temperature and 
E is the energy defined by Heff = — 'VE. Inserting this 
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equilibrium P{aj = 0,T) into Eq. (6), one finds that 



D 



1 



(7) 



this is the well-known dissipation-fluctuation relation. 
We now postulate that the fluctuating field is indepen- 
dent of the spin torque. This hypothesis is consistent with 
our notion that the spin torque is a deterministic action 
so that the spin torque does not alter the randomness 
induced by the thermal fluctuation. With this identifl- 
cation, the stochastic LL equation, Eq. (3), and its cor- 
responding Fokker-Planck equation, Eq. (6), completely 
determine the dynamics of the magnetization vector at 
finite temperature T. 



III. STATIONARY SOLUTION AND 
EFFECTIVE TEMPERATURES 

Before we numerically solve the above stochastic LL 
equation for a number of interesting cases, we should first 
look for a stationary solution in Eq. (6), i.e., P = PqCM.) 
is independent of time. Without the spin torque, this 
solution is known as the equilibrium Boltzmann distri- 
bution function mentioned earlier. With the spin torque, 
the system is no more in an equilibrium state because 
the system is subject to the spin torque and thus it is 
not a closed system. For an open system, the law of 
thermal dynamics does not require the minimum free en- 
ergy and the concept of thermal equilibrium breaks down. 
Nevertheless, it is still meaningful to obtain a stationary 
solution Ps where Fokker-Planck probability density is 
time-independent dPg/dt = 0. Thus, 



V ■ { [-7M X H,/^ + r, - X (M X H,//) 

-1-72(1 + a^)Dm X (M X V)] Ps} = 



(8) 



Unfortunately, the above eigenstate problem for an ar- 
bitrary field He// is generally difficult to solve. One 
can immediately verify that the Boltzmann probability 
density Pq oc exp(— P/fcsT) is no more a solution of 
the above equation. To make progress, we need to con- 
sider a special case below. First, we recall that the Neel- 
Brown formula for the thermal agitation is in fact most 
useful where the energy barrier constructed by He// is 
much higher than the thermal energy ksT . In this limit, 
the probability density will be very small if the direc- 
tion of the magnetization vector is away from the en- 
ergy minimum. Here we should also consider this case. 
We now tentatively seek a solution of P, in the form of 
Ps oc exp{—E/kBT*) where we have introduced an effec- 
tive temperature T* that will be determined next. By 
placing this P, into Eq. (8) and by utilizing the identity 
V • (M X He//P^) oc V • (M X VPs) = 0, we find 



V • { [a - 1) M X (M X He//) 
+ajM X (M X Mp)]Ps} ^ 



(9) 



Clearly, the above equation does not necessarily have a 
solution for an arbitrary effective field. However, as we 
point out earlier, we have limited ourselves to a high 
barrier case such that the magnetization vector at the 
stationary condition is nearly at the direction of Mp. 
For concreteness, let us consider a most experimentally 
relevant geometry where Mp ~ ex and 



H 



K 



77, 



Mx)ex - AuM^e, (10) 



where H^xt is the external field which is applied at x- 
direction, Hk is the anisotropy field, and —AirMzez is 
the demagnetization field perpendicular to the plane of 
the film. In this case, the energy minimum are at Mx = 



Mz = 0. We simply keep the first order in 
and set Mx = . I M'} 
to the first order, we find 



V ~ 
My and Mz, 



Ml 



Ml = Ms up 



M X (M X He//) w My{Hext + HK)ey 
+Mz{Hext + Hk + 47rA/,)ez 



(11) 



and M x (M x Mp) « Myey + M^Bz- By placing these ex- 
pressions into Eq. (9) and carrying out the divergence V 
for My and M^ components we find (one noticed that the 
divergence operator on P, produces higher orders since 
Eq. (11) is already the first order in My and M^, thus 
it is consistent with our approximation by neglecting the 
terms of divergence on P^ ) , 



a(27rMs + 



aj = 



or 



r* = T 1 - 



aj 



(12) 



(13) 



where we defined a critical spin torque Oc = a{H(,xt + 
Hk + ^ttMs). Coincidentally, this critical spin torque is 
precisely the minimum spin torque required to switch the 
magnetization at zero temperature [l£j. 

We should point out that the concept of the effective 
temperature introduced here should be understood in 
terms of the stationary solution of the probability density. 
The thermally averaged dynamical variable, for exam- 
ple, the magnetization vector < M >= J PgM. sin 6 dddcj) 
would behave as if the temperature of the system is 
T* . However, the magnetic temperature which is defined 
through the thermal fluctuation remains to be T. An al- 
ternative understanding of this effective temperature is to 
rewrite the stationary solution by Pq oc exp{—E* /kBT) 
where E* is an effective activation energy 



E* ^Eil 



aj_ 

0,0 



(14) 



In other word, we can state that the spin torque alters the 
magnetic energy and thus there will be an effective energy 
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1/KJ(mev)"' 

FIG. 1: Relaxation time r as a function of the inverse of the 
temperature for several values of the spin-transfer torques. 
The external field of H^xt = 440 (Oe) is applied along —x 
direction. 

barrier associated with the spin current. Therefore, it is 
equivalent to think of the effect of the spin torque via the 
modification of the temperature or of the energy barrier. 

To conclude this section, we have found a stationary 
solution Ps of the stochastic LL equation. Since the life 
time or the relaxation time r is inversely proportional to 
the probability density Ps, we can write the generalized 
Brown-Neel formula below 

r =/oexp(^ — j (15) 

where /o is an attempt frequency and Eh is the energy 
barrier. 



IV. COMPARISON BETWEEN NUMERICAL 
AND ANALYTICAL RESULTS 

The stationary solution, Ps oc exp[— i?(M)/fcBT*], is 
based on the assumption of high energy barrier assump- 
tion. In general cases, one should start the calculation of 
magnetization dynamics from our generalized LL equa- 
tion, Eq. (3). Since the stationary solution is simple and 
easy to use in analyzing experimental data, it would be 
necessary to establish the range of validity of Eq. (15) 
for various interesting experimental situations. Once its 
validity is established, we expect that our stationary solu- 
tion would be serving as a starting point to understand 
various thermal agitation phenomena in the present of 
the current. In this Section, we numerically integrate 
Eq. (3) and compare the result with Ps- 

The standard white-noise spectrum, Eq. (4) where D 
is given by Eq (7), is used for the modeling of the tem- 
perature dependence of random fields. The calculation 
procedure is same as that for the standard LLG equa- 
tion, except a spin torque is added to the equation of 
the motion. A magnetic layer, whose lateral size is 
64nmx64nm and whose thickness is 2.5nm, is treated 



as a single macrospin. The in-plane uniaxial anisotropy 
field Hk is 500(Oe) and the saturation magnetization is 
AttMs — 12,000 (Oe). These parameters are reasonably 
consistent with the experiments performed by the Cornell 
group The Gilbert damping constant was taken as 
a = 0.03 throughout the modeling. The magnetization 
of the free layer is initially saturated at -I- a; direction. At 
t = 0, we apply a magnetic field at —x direction whose 
magnitude is close to but less than the anisotropy field 
Hk- At the same time, aj is also applied to the system. 
Equation (3) was numerically integrated in time using 
the stochastic Heun method with a 0.3ps time step. A 
smaller time step has been tested and it yields nearly 
identical results in all the cases presented in the paper. 

With above specified parameters and procedure, we 
first determine the probability P{t) that the magnetic 
layer has been reversed within the waiting time t. By 
performing up to 5 x 10^ independent runs for each set of 
parameters (each run starts at i = and ends at the time 
that the magnetization has either been just switched or 
ends at the time up to t = 5/is, whichever is smaller) and 
then by recording the number of them that the magneti- 
zation is switched at time interval (i, t+At), we obtained 
a simulated P{t) that is fitted by a simple exponential 
function, i.e., P{t) = 1 — exp(— f/r) where r is the fitting 
parameter for the relaxation time. We have found that 
the fit works remarkably well for any values of a j we have 
considered. In Fig. 1, we show the fitted relaxation time 
T as a function of the temperature for a fixed applied 
field Hext — 440Oe. Two features are immediately seen: 
first the data fall on a straight line for any fixed aj; this 
indicates the thermally assisted reversal can be described 
by an activation process, i.e., Inr cx l/ksT. Therefore, 
it is meaningful to introduce an effective activation en- 
ergy, see Eq. (14), in accordance with the Neel-Brown 
law of thermal agitation. The second conclusion is that 
the slope, or the activation energy depends on STT: the 
positive aj favors a lower energy barrier. All these fea- 
tures are well described by Eqs. (13)-(15). In Fig. 2, the 
effective activation energies are shown to be linearly de- 
pendent on the current and they vanish at almost the 
same point ac for different external fields [note that Oc 
is weakly dependent on the external field, see the defini- 
tion of Oc after Eq. (13)]. In the insert of Fig. 2, we have 
shown the activation energy as a function of the magnetic 
field for several different STT. The activation energy can 
be fitted by 



where Hs is the switching field at zero temperature, Eq 
is the energy barrier at zero magnetic field field, and (3 
is a constant, which has been argued to be 1.5 or 2. The 
exponent (3 — 2 ior the external field applied parallel 
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Spin transfer torque (1 0^ Oe) 

FIG. 2: Activation energy as a function of aj at three 
different external fields. Inset: El vs 1 — Hsxt/Ha for several 
different values of STT. 



to the easy axis. These simulated results confirm our 
analytical result, Eq. (14). 



V. COMPARISON WITH EXPERIMENTAL 
DATA 

A number of experiments on spin torque induced ther- 
mal agitation had been carried out. It would be inter- 
esting to see whether our prediction, Eq. (15), agrees 
with these existing data. The phenomenon that we want 
to compare first is so-called "telegraph noises" or dwell 
times. Experimentally, one simultaneously applies an ex- 
ternal magnetic field and a spin current to a spin valve 
structure so that the magnetization direction of the free 
layer is fluctuating from one direction to another due to 
thermal agitation The dwell time rp (tap) is 

defined as an average time the magnetization of the free 
layer is parallel (antiparallel) to that of the fixed layer. In 
general, rp ^ tap- However, by adjusting the magnetic 
field or the spin current, one is able to obtain an equal 
dwell time for parallel and antiparallel states, rp — tap- 
From Eq. (15) for the parallel and the antiparallel states, 
the condition of the equal dwell time is 

(17) 

where I^^ and are the critical currents for the mag- 
netization switching from antiparallel to parallel align- 
ments and vice versa; their magnitudes are not necessary 
the same, i.e., I^^ ^ —I^ in a typical experimental ge- 
ometry [13, To compare our prediction, Eq. (17), 
with experiments, we plot the H-I phase diagram of equal 
dwell time rp = tap in Fig. 3. 

It is noted that we have shifted the external field by 
= — 22(mT) to take into account of the magnetic 
coupling between the free and fixed layers. The cou- 



20- 




-100^ • ^ ■ ^ ■ ^ ■ 1 

-0.6 -0.3 0.0 0.3 0.6 

Current I (mA) 



FIG. 3: H-I Phase boundary of equal dwell times < tap >~< 
TP >. The coupling field is Hah = —22mT, the critical cur- 
rents are I^^ = —1.25mA for the transition from AP to P 
alignments and = O.A25mA from P to AP alignments. 
Line: Eq. (17) except He^t being replaced by H^xt + Hsh- 
Square: experimental data [l^. 




Current I (mA) 

FIG. 4: Ratio of the relaxation times ln(/orp)/ ln(/oTAp). 
The parameters are same as those in Fig. 4. 

pling may come from either the exchange or dipolar cou- 
plings. Since we assume that the magnetization of the 
fixed layer is held at the direction of e^,, the free layer 
receives an effective coupling field that will be added to 
Hext in Eq. (17). Evidently, the agreement between our 
results and the experimental data is excellent 0, 0| . 

Next we compare the ratio of the dwell times of an- 
tiparallel and parallel states for a fixed magnetic field. 
Again, from Eq. (15), we have 



ln(/orp) 
ln(/oTAp) 




In Fig. 4, we show the ratio of the dwell times for a 
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T T=300 K 



0.001 0.002 0.007 0.018 0.050 

C = daj/dt (Oe/ns) 

FIG. 5: Dependence of ac(T,C) on the sweeping rate C at 
finite temperature. Inset: ac{T,C) and Ua as a function of 
temperature for the sweeping rate of STT at O.OlOe/ns. 



fixed magnetic field as a function of the spin torque by us- 
ing the same set of experimental parameters as in Fig. 3. 
The results are consistent with experimental data (how- 
ever, the data points in Ref . [13, are rather scattered 
so that I do not include those data in the figure). 

Up till now, we have studied the thermal activation 
by abruptly introducing an external field and a STT at 
t=0. In experiments, there are ramping times, e.g., STT 
is gradually increasing at rate of, say 10~^Oe//is ^3 
and care must be taken when one compares our theory, 
Eq. (15), with experiments. In the current ramping pe- 
riod, a J is not a constant and thus the activation energy, 
Eq. (16), varies with time. In this case, one should utilize 
the differential form of the switching probability instead 
oiP{t) = 1 -exp(-i/r). 



dP{t) _ dt 



(19) 



Numerically, the finite ramping rate can be rather eas- 
ily handled. In determining ac{T,C) at finite tempera- 
ture, we ramp aj with a fixed rate. At a certain time, the 
magnetization vector switches and we record the value of 
aj. By repeating the above procedure 800 times, we are 
able to establish the switching aj histograms from which 
the mean switching ac{T,C) and its standard deviation 
Ua are obtained. At temperature between lOO-ftT and 
SOOi^T, we calculated the mean and standard deviation of 
the distributions at a sweeping rate between O.OOlOe/ns 
and Q.lOe/ns. As expected for a thermally activated 
process, ac(T, C) increases with decreasing temperature 
and with increasing sweeping rate. Figure 5 shows tem- 
perature and sweeping rate dependence of ac{T,C) and 

CTa- 

A logarithmic dependence of adT, C) on the sweeping 
rate has been found, which is in a good agreement to 
Eq.(20). Moreover, the inset of Fig. 5 describes the tem- 
perature dependence of ac(T, C) and Oa- We have ver- 
ified that the temperature dependence of CTq is a linear 
relationship and adT^C) monotonically decreases with 
increasing T . 

Myers et al have discovered that the thermal ac- 
tivation driven by spin torque is qualitatively different 
from that driven by the magnetic field. They have sug- 
gested an activation energy whose form is similar to ours, 
except that they have postulated an arbitrary exponent, 
i.e., 5W oc (1 — aj/ttc)^. Although they have stated that 
^ might be 1.5, most of experimental data shown in the 
paper can be used to determine the value of ^. One set of 
data. Fig. 1(d) of Ref. [11], shows that ac{T,C) linearly 
increases with InC as predicted by our Eq. (15). If one 
uses different scaling relation, e.g., SW oc (1 — aj /a^Y'^ , 
one would obtain ac(T,C) oc (InC)^^^ that would dis- 
agree with the experimental data. Therefore, the exist- 
ing data support the linear scaling between the activation 
energy and the spin torque. 



dai 



By assuming a linear ramping of STT, i.e., C — is a 
constant and by placing Eqs. (15) and (16) into Eq. (19), 
we integrate Eq. (19) from t — to t — to and find the 
average switching STT ac{T, C) = aj{ta) 



ac{T, C) = flc 

where -/oln[l- 
ing STT is found as 



EbC 



(20) 



P(to)]- The variance of the switch- 



' Eh 



(21) 



These relations, Eqs. (20) and (21), are consistent with 
earlier studies on the similar energy barrier formalism for 
a very different physical system |18| . 



VI. CONCLUSIONS 

In summary, we have extended the law of thermal ag- 
itation to include the spin transfer torque driven by the 
spin polarized current in magnetic multilayers. Although 
the concept of the energy barrier or the temperature in 
the Neel-Brown formula breaks down in the presence of 
the spin transfer torque, we are able to re-establish the 
Neel-Brown formula by properly introducing an activa- 
tion energy or an effective temperature to replace the 
true energy barrier or true lattice temperature. Our for- 
malism is further supported by numerical solutions and 
is in agreement with experimental results. 

This work is supported by NSF (ECS-0223568). 
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